Direct Josephson coupling between superconducting flux qubits 
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We have demonstrated strong antiferromagnetic coupling between two three- j unction flux qubits 
based on a shared Josephson junction, and therefore not limited by the small inductances of the qubit 
loops. The coupling sign and magnitude were measured by coupling the system to a high-quality 
superconducting tank circuit. Design modifications allowing to continuously tune the coupling 
strength and/or make the coupling ferromagnetic are discussed. 

PACS numbers: 74.50.+r, 85. 25. Am, 85.25.Cp 



Quantum superposition of macroscopic states was con- 
clusively demonstrated in superconducting Josephson 
structures in 2000. 1 Such structures are natural candi- 
dates for the role of qubits (quantum bits), the con- 
stituent elements of quantum computers. Successful op- 
eration of a quantum computer would be the ultimate 
confirmation of the validity of quantum mechanics on the 
macroscopic scale, which makes the task of controllably 
linking a significant number of qubits together more than 
just an advance in technology. 

The coupling energy J must be comparable to the split- 
tings between the two lowest eigenstates of individual 
qubits. On the other hand, the coupling must not excite 
the qubits to higher levels, or significantly increase the 
qubits' interaction with undesirable degrees of freedom, 
leading to decoherence and dissipation. Finally, J should 
be either variable by design or, even better, tunable dur- 
ing the system's operation. 

In this letter we demonstrate the coupling of two three- 
Josephson-j unction (3 J J) flux qubits, making progress to- 
wards meeting these requirements, and discuss the ways 
of its further improvement. The 3JJ qubit consists of 
a superconducting loop with small inductance L inter- 
rupted by three Josephson junctions. The two differ- 
ent directions of persistent current in the loop form the 
qubit 's basis states. 2 The 3JJ design enables classical 
bist ability even for L — > 0, resulting in a weak coupling 
to environmental magnetic flux noises. As a result, quan- 
tum behaviour with long decoherence time was observed 
in this type of qubit by several groups. 3 " 5 

However, their small L makes it difficult to couple 
3JJ qubits inductively; generally, J is smaller than the 
single-qubit level splitting. We therefore implement the 
proposals 6 ' 7 to directly link two qubits through a shared 
junction (Fig. 1). The resulting coupling not only is 
strong, but can also be varied independently of other de- 
sign parameters by choosing the shared junction's size. 




FIG. 1: Schematics of a 7JJ device: two 3JJ qubits with direct 
Josephson coupling and reduced flux bias f a -b- 



To calculate J, we neglect the inductances so that 
the potential term in the Hamiltonian contains only 
the Josephson energy Uj = — J2^=o Ej cos 0j •> an d use 
flux quantization, 0i + 02 + 03 + 0o = 27r(^+/ a ) and 
04 + 05 + 06 - 0o = 27r(!+/ 6 ), to eliminate 2;5 (f a; b = 
^a-b/^o — \ ^ S tne reduced flux bias). In the simplest 
case \fa-b = 0; £a ; 3 ; 4 ; 6 = E\ E 2 -$ = aE (| < a < 1); 
Eq ^> E] there are two different pairs of potential min- 
ima: "ferro-" and "antiferromagnetic" (with parallel and 
antiparallel loop currents respectively), 
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where 7 p = (2e/H)E^/l — l/4a 2 is the persistent current 
in a free 3JJ qubit. 2 Inserting these into [7j, one finds 
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that the AF states have the lower energy by 
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so that the effective mutual inductance h 2 /4e 2 Eo is 
just the standard Josephson inductance of the coupling 
junction. 8 For an explanation, note that flipping the signs 
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yields an AF configuration with 



0, and 



with the same energy as the FM minimum. Such a state 
of course is non-stationary (since charge must be con- 
served), and adjustment of the phases will lower [7j, with 
the stationary AF state (2) realizing the global minimum. 
More intuitively: the nonzero c/)q F reduces the effective 
frustration in the individual qubits, which is maximal 
for f a .b = [cf. above (1)]; the attendant reduction in 
qubit energy overcomes (by a factor two) the increase in 
Josephson energy in the coupling junction itself. 

Thus the direct Josephson coupling of 3JJ qubits has 
the same sign as their inductive coupling (the latter cor- 
responding to the natural north-south alignment of their 
magnetic moments), 9 but is not restricted by the geomet- 
ric inductances. For Eq — > oo, J in (3) disappears as it 
should, since then we have two qubits sharing a com- 
mon leg without a junction — equivalent to two adjacent 
qubits if kinetic inductance is neglected, with only induc- 
tive coupling. 10 

We determined J using an impedance measurement 
technique, applied previously to 3JJ qubits 4 ' 11 and ex- 
tended to multiple qubits in Ref. 12. The qubits 
are placed inside a tank circuit with known self- 
inductance L T and quality factor Q T , driven by a dc bias 
plus a small ac current at its resonance frequency cut 
(Fig. 2). The tank's voltage-current phase angle is 
given by tan# = — (Qt/Lt)x'- Here, \' is the qubits' 
contribution to the tank susceptibility, 13 readily related 
to the curvature of their energy bands: for qubit-tank 
mutual inductances M a x = M b T = M, one has 14 
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where <£ x denotes a symmetric change of flux bias in both 
qubit loops. At temperature T = 0, E tot is the qubits' 
ground-state energy; at finite T, the derivative simply 
becomes a Boltzmann average over the levels. The band 
curvature is large near anticrossings, so that tan0(/ a , fa) 
contains important information about the level structure. 

One obtains J from such a plot as follows. The stan- 
dard four-state Hamiltonian for two coupled qubits is 



H 



e b a z b - A b <j x b + Ja z a z b , (5) 



where o x , cr z are Pauli matrices, Aj is the tunnel- 
ing amplitude, and ej = I P j$ofj is the energy bias 
(j = a, b). For low T and small Aj, the location of 
the peaks in |tan0| (due to anticrossings) follows sim- 
ply from the classical stability diagram. For instance, 
the |t|) <-> |tt) transition (|TI)- o~ z = -a b = 1 etc.) oc- 
curs at —e a + e b — J = —e a — e b + J e b = J. Therefore, 
the peak-to-peak distance equals 2 J. 



FIG. 2: Micrograph of sample 1 (see the table). The 3JJ 
qubits a and b are coupled through a shared Josephson junc- 
tion, visible in the centre, and can be flux-biased indepen- 
dently by dc lines ibi and ib2- The pancake coil is part of an 
LC circuit reading out the qubits' magnetic susceptibility. 



For our samples, we first fabricate niobium (Nb) pan- 
cake coils and dc flux-bias lines on 4-inch oxidized silicon 
wafers, and then the qubits inside the coils by aluminium 
(Al) shadow evaporation on 12x12 mm 2 chips. 

The Nb process starts with sputtering and dry etching 
of the 200 nm thick coil windings with 1 /am width, 1 /im 
line spacing, and typically 30 turns. The patterning uses 
e-beam lithography and a CF4 RIE process. Then a sil- 
icon oxide isolation layer and the second 300 nm thick 
Nb film are deposited for the central coil electrode and 
the 2 /im wide dc lines; photolithography is used for all 
required resist masks of these layers. Finally, 400 nm 
silicon oxide is deposited for protection and isolation. 

The Al process uses e-beam lithography to prepare the 
double layer resist mask for the qubits with a 150 nm 
linewidth. The two Al layers are deposited in situ by 
e-beam evaporation with different angles of incidence at 
a rate of 1.8 nm/s. The surface of the first Al film is 
oxidized with pure oxygen at a pressure of 10 -2 mbar. 
The qubits are completed after the final lift-off. 

Results are shown in Fig. 3(a). As explained below (5), 
one has J « I pa ®o8fa = Ipb^odfa- One can find I pa;b in 
two different ways , which ag ree to within ~20%. First, 
we used I p « I c ^/l — l/4a 2 [cf. below (2)]. Here, I c is 
the critical current of a junction fabricated on the same 
chip and with the same area of 650 x 150 nm 2 as junctions 
1/3/4/6, enclosed in a superconducting loop and mea- 
sured by the conventional rf- SQUID technique; a w 0.75 
by design. The second way is to fit the shape of the peaks 
in taxi 6(f a , fa) [Fig- 3(b)], using the spectrum of (5) to 
evaluate (4), which yields Aj and ipj. 11 ' 12 ' 15 The re- 
quired tank parameters were extracted from its resonance 
characteristic, and the mutual inductances from the flux 
periodicity; for sample 1, Lt = 136 nH, Qt = 664, 
u t /2tt = 19.925 MHz, and M aT = M bT 66.5 pH. 
The |t|) <-> ||t) anticrossing does not show up in the 
figure because there is no net flux change, hence no con- 
tribution to the qubit susceptibility; one can also say that 
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FIG. 3: Plot of — tan 0(/ a ,/b) for two coupled 3JJ qubits 
(sample 1 from the table), (a): Measurements at a nominal 
mixing- chamber temperature of 20 mK (actual data taken 
between the dashed lines) . (b) : Theoretical fit for T = 70 mK 
(see the text for the different T's). The coupling strength can 
be estimated from the separation between the peaks. 



TABLE I: Coupling-junction areas So, tunneling ampli- 
tudes Aj, persistent currents I p j, coupling energies J, and 
peak locations Sf for the measured samples. 



sample 
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A a 


A, 
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Ipb J <5fexp 5fth 


No. 
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mK 


mK 


nA 


nA K 


1 


0.30 


80 


90 


120 


110 0.7 0.037-0.041 0.0360 


2 


0.15 


30 


30 


150 


120 1.2 0.068 0.0675 



the level curvature peaks perpendicular to the symmetric 
direction stipulated in (4). 

The results of the fit are summarized in the table for 
two measured two-qubit samples, with different sizes of 
junction 0. Note how, say, A a < A5 for sample 1 makes 
the a-anticrossing sharper, resulting in a deeper colour 
for the corresponding peak (vertical bands in Fig. 3). 
Since Eo/E = 3.1 and 1.5, respectively, the pertur- 
bative analysis leading to (3) does not apply quanti- 
tatively (the latter would have required large coupling 
junctions, which proved difficult to make with sufficient 
homogeneity). Instead, a theoretical prediction is made 
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FIG. 4: As in Fig. 3, but for a mixing-chamber temperature 
of 200 mK (a) and an effective T = 300 mK (b). 



for 5f(Eo/E,a), by calculating the classical stability di- 
agram directly from f7j(0o, </>i, 03, </>4, (j>s)- Incidentally, 
this has the advantage that the critical-current density 
drops out of the comparison (last two columns in the ta- 
ble), which therefore shows greater accuracy than we can 
claim for J itself. 

Data taken at a higher T support the effective Hamil- 
tonian (5) for our 7JJ system beyond the ground state. 
Namely, for, say, f a < 0.035, f b w 0.042, an |U) ^ | TT> 
anticrossing persists between excited states of sample 1. 
At finite T, it should contribute in (4), with rapidly de- 
creasing Boltzmann weight as f a is reduced. This is pre- 
cisely what is seen in Fig. 4(a); the fit in Fig. 4(b) shows 
detailed agreement with the theory. In both Figs. 3 and 4, 
the discrepancy between effective and mixing-chamber 
temperatures is well within the range expected due to 
heating through external leads etc.; we observed no sig- 
nificant deviations from an equilibrium distribution. 

The remarkable J ~ 1 K significantly exceeds both the 
tunnel splitting and the inductive coupling (estimated to 
be ~20 mK). It can be flux-tuned by using a standard 
compound junction (dc-SQUID) for the coupling. In- 
stead, one can also apply a bias current 1^ through junc- 
tion 0. 16 The corresponding generalization of (3) reads 
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J = hll/2eyjr^-T* + O[(E/E ) 2 }. Thus, J can only 
be increased, albeit significantly. Hence, this mecha- 
nism does not allow, e.g., changing the coupling sign 
and tunable decoupling of qubits. These are desirable 
for most quantum algorithms, but existing proposals for 
flux qubits rely on, and therefore are limited by, mu- 
tual inductances. 17 A bias line will introduce some noise. 
For reference, we give the coupling linewidth due to low- 
frequency fluctuations in 1^ with spectral density Sb(u): 
AJ = M p 4 /2 5b (0)/4e 2 (/ c 2 o -/ b 2 ) 3 . 




FIG. 5: Two 3JJ qubits with ferromagnetic Josephson cou- 
pling, (a): Twisted design, (b): Overlapping design. 



Other variations are presented in Fig. 5. In Fig. 5(a), 
the relative twist between the qubit loops interchanges 
the role of the AF and FM configurations, so that the 
latter have the lowest energy. In particular, the strength 
of the direct FM Josephson coupling can overcome any 
residual AF inductive interaction. Junction can pre- 
sumably be fabricated between the crossing lines in the 
centre. In Fig. 5(b), one qubit loop is 01230 and the 
other is 04560; by choosing 1:2 area ratios, both qubits 
can be brought close to degeneracy with a homogeneous 
field, for = |<E>£ = ^<f> . One obtains FM coupling 
without a twisted layout, but with strongly asymmetric 
qubits. The two discussed modifications can be com- 
bined: by current-biasing the junctions in Fig. 5, one 
obtains tunable FM coupling. 

In conclusion, we have demonstrated direct antiferro- 
magnetic Josephson coupling between two individually 
controllable three-junction flux qubits. The coupling 
strength can be on the order of a Kelvin. We also pro- 
posed design modifications allowing tunable and/or fer- 
romagnetic coupling. Future experimental work should 
also consider linear qubit arrays, to which the design of 
Fig. 1 is readily generalized. 7 
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